Universal dynamical decoherence control of noisy single- and 

multi-qubit systems 

Goren Gordon*, Noam Erez^, Gershon Kurizki* 
Department of Chemical Physics, Weizmann Institute of Science, 76100 Rehovot, Israel 

Abstract 

In this article we develop, step by step, the framework for universal dynamical control of two- level 
systems (TLS) or qubits experiencing amplitude- or phase-noise (AN or PN) due to coupling to a 
thermal bath. A comprehensive arsenal of modulation schemes is introduced and applied to either 
AN or PN, resulting in completely analogous formulae for the decoherence rates, thus underscoring 
the unified nature of this universal formalism. We then address the extension of this formalism to 
multipartite decoherence control, where symmetries are exploited to overcome decoherence. 
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I. INTRODUCTION 



In-depth study of the mechanisms of decoherence and disentanglement and their preven- 
tion in bipartite or multipartite open systems is an essential prerequisite for applications 
involving quantum information processing or communications [lj. The present article is 
aimed at furthering our understanding of these formidable issues, which is scanty at best. It 
is based on recent progress by our group, as well as others, towards a unified approach to the 
dynamical control of decoherence and disentanglement. This unified approach culminates 
in universal formulae allowing design of the required control fields. 

The topic of multipartite decoherence has been well-investigated in two limits. One of 
these is relaxation toward steady-state of one-body coherence of spins, atoms, excitons, 
quantum dots, etc. that are in contact with a much larger reservoir. The other is the 
collective decoherence of a small (localized) two-body or many-body system, which typically 
occurs more rapidly than one-body decoherence j^, 3]. 

By contrast, more general problems of decay of non-local mutual entanglement of two 
or more small systems are less well understood. This decoherence process may occur on a 
time scale much shorter than the time for either body to undergo local decoherence, but 
much larger than the time each takes to become disentangled from its environment. The 
disentanglement of individual particles from their environment is dynamically controlled by 
interactions on non-Markovian time-scales, as discussed below 0]. Their disentanglement 
from each other, however, may be purely Markovian , in which case the present 

non-Markovian approach to dynamical control/prevention is insufficient. 

A. Dynamical control of decay and decoherence on non-Markovian time scales 

Quantum-state decay to a continuum or changes in its population via coupling to a ther- 
mal bath is known as amplitude noise (AN). It characterizes decoherence processes in many 
quantum systems, e.g., spontaneous emission of photons by excited atoms |2|], vibrational 
and collisional relaxation of trapped ions 8] and the relaxation of current-biased Josephson 
junctions [9|. Another source of decoherence in the same systems is proper dephasing or 
phase noise (PN) [3| , which does not affect the populations of quantum states, but random- 
izes their energies or phases. 
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A thoroughly studied approach to suppression of decoherence is the "dynamical decern- 
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pling" of the system from the bath jlO . 
"bang-bang" (BB) pulses have been proposed for stroboscopic suppression of proper dephas- 
ing: 7r-phase flips of the coupling via strong and sufficiently fast resonant pulses applied to 
the system [3, 14, 15]. The identification of a decoherence-free subspace (DFS), wherein 
symmetrically degenerate states are decoupled from the bath, constitutes a complementary 

3- 

Our group has purported to substantially expand the arsenal of decay and decoherence 
control. We have presented a universal form of the decay rate of unstable states into any 
reservoir (continuum), dynamically modified by perturbations with arbitrary time depen- 
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dence, focusing on non-Markovian time-scales 

been obtained by us for the dynamically modified rate of proper dephasing 
Our unified, optimized approach reduces to the BB method in the particular case of proper 
dephasing or decay via coupling to spectrally symmetric (e.g., Lorentzian or Gaussian) noise 
baths with limited spectral width (see below). The type of phase modulation advocated for 
the suppression of coupling to asymmetric baths (e.g., phonon or photon baths with fre- 



quency cutoff 



28]) is, however, drastically different from the BB method. Other situations 



to which our approach applies, but not the BB method, include amplitude modulation of the 
coupling to the continuum, as in the case of decay from quasibound states of a periodically 
tilted washboard potential 24j: such modulation has been experimentally shown 29J to 
give rise to either slowdown of the decay (Zeno-like behavior) or its speedup (anti-Zeno-like 
behavior), depending on the modulation rate. 

The theory has been generalized by us to finite temperatures and to qubits driven by an 
arbit rary time-dependent field, which may cause the failure of the rotating-wave approxima- 
tion [25(] . It has also been extended to the analysis of multi-level systems, where quantum 
interference between the levels may either inhibit or accelerate the decay 301 ] . 

Our general approach 24j to dynamical control of states coupled to an arbitrary "bath" 
or continuum has reaffirmed the intuitive anticipation that, in order to suppress their decay, 
we must modulate the system-bath coupling at a rate exceeding the spectral interval over 
which the coupling is significant. Yet our analysis can serve as a general recipe for optimized 
design of the modulation aimed at an effective use of the fields for decay and decoherence 
suppression or enhancement. The latter is useful for the control of chemical reactions 31]. 
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B. Control of symmetry-breaking multipartite decoherence 



Symmetry is a powerful means of protecting entangled quantum states against deco- 



herence, since it al 



subsystem 



ows the existence of a decoherence-free subspace or a decoherence-free 
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34J. In multipartite 



systems, this requires that all particles be perturbed by the same environment. In keeping 
with this requirement, quantum communication protocols based on entangled two-photon 
states have been studied under collective depolarization conditions, namely, identical ran- 



dom fluctuations of the polarization for both photons 



351 ] . Entangled qubits that reside at 



the same site or at equivalent sites of the system, e.g. atoms in optical lattices, have likewise 
been assumed to undergo identical decoherence [4|. 

Locally- decohering entangled states of two or more particles, such that each particle 
travels along a different channel or is stored at a different site in the system, may break the 
state symmetry. A possible consequence of this symmetry breaking is the abrupt "death" 
of the entanglement Such systems, composed of particles undergoing individual 

or "local" decoherence, do not possess a natural DFS and thus present more challenging 
problems insofar as decoherence effects are concerned [36 ]. 

Our group has recently addressed these challenges by developing a generalized treatment 
of multipartite entangled states (MES) decaying into zero-temperature baths and subject to 
arbitrary external perturbations whose role is to provide dynamical protection from decay 
and decoherence [30, [3T|] - Our treatment applies to any difference between the couplings 
of individual particles to the baths. It does not assume the perturbations to be strobo- 
scope, i.e. strong or fast enough, but rather to act concurrently with the particle-bath 
interactions. Our main results are to show that by applying local (selective) perturbations 
to multilevel particles, i.e. by addressing each level and each particle individually, one can 
create a decoherence-free system of many entangled qubits. Alternatively, one may reduce 
the problem of locally decohering MES to that of a single decohering particle, whose dy- 
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261 ] . On the other hand, 



namical control has been thoroughly investigated 
the combined effect of dephasing and relaxation (phase and amplitude noises) on MES and 
its control, which constitute a much more formidable problem, have not yet been studied by 
us. 
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C. Outline 



In this article we develop, step by step, the framework for universal dynamical control 
by modulating fields of two-level systems or qubits, aimed at suppressing or preventing 
their noise, decoherence or relaxation in the presence of a thermal bath. To this end, a 
comprehensive treatment is developed in Sec. [IT] hi a more complete and transparent fashion 
than its brief sketch in Ref. [251]. Its crux is the derivation of a more general master equation 
(ME) than in previous treatments of a multilevel, multipartite system, weakly coupled to an 
arbitrary bath and subject to arbitrary tem por al driving or modulation. The present ME, 



derived by the Nakajima-Zwanzig technique 38|, |39j], is more general than the ones obtained 



previously in that it does not invoke the rotating wave approximation and therefore applies 
at arbitrarily short times or for arbitrarily fast modulations. 

Remarkably, when our general ME is applied to either AN or PN in Sec. 1111} the resulting 
dynamically-controlled relaxation or decoherence rates obey analogous formulae provided 
the the corresponding density-matrix (generalized Bloch) equations are written in the ap- 
propriate basis. This underscores the universality of our treatment. The choice of an ap- 
propriate time-dependent basis allows here to simplify the AN treatment of Ref . 25] . More 
importantly, it allows us to present a PN treatment that does not describe noise phenomeno- 
logically as in Ref. 25j, but rather dynamically, starting from the ubiquitous spin-boson 
Hamiltonian. 

We then discuss in Sec. IIV} more comprehensively than in previous treatments, the pos- 
sible modulation arsenal for either AN or PN control. The present formalism is appli- 
cable in a natural and straightforward manner to multipartite and/or multilevel systems 



37| | . It allows us to focus in Sec. |V] on the ability of symmetries to overcome multipartite 



decoherence 



18 



19 



20, |2l|, |34|. Our conclusions are presented in Sec. ED 
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II. MASTER EQUATION (ME) FOR DYNAMICALLY CONTROLLED SYS- 
TEMS COUPLED TO THERMAL BATHS 



A. Derivation of the reduced density matrix ME by the Nakajima-Zwanzig 
method 

We shall consider the most general unitary evolution of a system coupled to a thermal 
reservoir, governed by the Liouville operator equation (we shall take H — 1 throughout the 
rest of the paper): 

p tot (t) = -i[H(t),p tot (t)] = -iL(t)p tot (t), (1) 

H(t) = H (t) + Hj(t), (2a) 
H (t) = H s (t)+H B , (2b) 

L(t) = L s (t)+L B + Lj(t) (3) 

where p to t is the density matrix of the system+reservoir and if 5, Hb, Hj are the Hamiltoni- 
ans of the system, bath and their interaction, respectively. As usual, L denotes the Liouville 
operator, which acts linearly on operators on our Hilbert space. We shall use the notation 
Ho(t) for the "unperturbed" Hamiltonian, assuming weak system-bath coupling, and Lo(t) 
for its associated Liouville operator. 

We seek a master equation for the reduced density matrix of the system alone p = 
TrsPtoti allowing for arbitrary time dependence of H(t) without resorting to the rotating 
wave approximation 

as 

. This can be acc omp lished by the Nakajima-Zwanzig S. v 
projection-operator technique (see also 4l|, |42j). Let us define ps = Z~ 1 e~ HB / kBT (T being 
the reservoir temperature, Z normalization to unit trace), the projection operator V(-) = 
Ttb(') ® Pb (satisfying V 2 = V), and the complementary (projection) operator Q = 1 — V. 
In terms of these definitions, Eq.flT]) is equivalent to: 

V Ptot (t) = -tVL(t)V Ptot (t)-iVL(t)Q Ptot (t) (4a) 
Qptotit) = -tQL{t)V Ptot {t)-iQL{t)Qp tot {t) (4b) 



Equation (14bl) is then formally integrated to give: 



Qfh*(t) = -i [ JC + (t,T)QL{T)V Ptot (T)dT + )C + {t,0)Q Ptot (0) 
Jo 



(5) 



where 



JC+(t,r) = T + e^ s ^ L(s)ds , (6) 

T + denoting time-ordering. 

If this expression is then plugged into Eq. (j4aj) we get a non-Markovian ME for Pp to t- 



V Ptot (t) = -iVL(t)V Ptot (t)- [ VL(t)lC + (t,T)QL(T)Vp tot (r)dr-tVLIC + (t,0)Qp tot (0) (7) 

Jo 



which yields a ME for p after tracing out the bath. 

Rather than apply a perturbative treatment directly to t 
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lis equation, it is useful to first 
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461 ]. In this form, the 



transform it to a "time-convolutionless" form (TCL) 
memory effect (the presence of Ptot( r ) i n the integrand) is transferred to the integration 
kernel and p(t) is taken out of the integral. Only then, is the perturbative expansion (in hi) 
applied. 
Formally, 



Ptot(r) = G-(t,r)p tot (t) 
where, writing T_ for anti-chronological ordering: 

£_(t,r) = T_e + ^ L(s)ds . 
Substituting this expression for p to t( T ) into Eq.Q, one obtains: 



(8) 



(9) 



Qptot(t) = - [ K + (t,T)iQh(T)Vg_(t,r)dT(V + Q)p to t(t)+}C + (t,0)Qp tot (0) (10a) 
Jo 

Collecting all Qptot terms on the left, we obtain: 



Ht) Qptot (*) = {!- F(t) } Vptot (t) + K + (t, 0) Q Ptot (0) 

T{t) =1+ I ]C + (t,T)iQL(T)VG-(t,T)dT = l + E(t) 
Jo 



(10b) 
(10c) 



Assuming T{t) can be inverted (which is expected to hold for short times in the weak 
coupling limit), and writing 0(t) = jF(t) -1 , one obtains the equation: 

Qptotit) = {&(t) - l}V Ptot (t) + 6(t)/C + (t,0)Qp tot (0). (lOd) 
Finally, plugging this expression for Qp to t into Eq. (|4aj) . the formal TCL ME is obtained: 

Vp tot (t) = -iVL(t)V Ptot (t)-iVh{Q(t) - l}V Ptot (t) -zPL0(t)/C + (t,O)Q Ptot (O). (11a) 

If the initial condition is such that Qptot(O) = 0, so that the last term vanishes (as will 
indeed be our assumption below), then all memory effects are contained in O(t). In what 
follows, we shall always assume that 

Ptot(0) = p s (0)®p B , (lib) 

so that this condition is fulfilled. 

The operators L5 and L B both commute with V and Q, and VL B = 0. This implies 
VLQ = VhjQ (and note that - 1 = Q(0 - 1)). With the notation {-) B = Ti B (-p B ), the 
ME for the reduced density matrix of the system can be written in the form: 

P (t) = -1 [L s + (L/) B ] pit) - E{t)p(t) (12a) 
E(t) = (iLi{Q{t)-l}) B (12b) 

B. Born approximation 



43, 



48J]. We begin 



At this point, it is expedient to follow the perturbational method of 
by noting that (cf. Eqs.ffT0a.ffT2bll): 



S(t) = -2(L 7 {S(t) + 0(S 2 )}) B ; E(t) = 0(Lj). (13) 

Therefore, to expand Eq. fll2bj) to second order in Lj, we need to evaluate /C + and G only to 
0th order! The operator /C+(t, r) can be factored: 

/C + (t, r) = V (t, T)T + e- iQ ti Vo ^ lL '^ Vo ^ ds (14a) 



S 



where 



V (t,T) = T +e - iS ^ Lo(s)ds . (14b) 
To zeroth order in Lj, fC + is just Vq. Furthermore, it is important to note that 

V Q = QU Q = QU (15a) 

where 

U (t,r) = T + e- i £ L °W da , (15b) 
and in Eq.f ll5ap we have used the fact that Uq commutes with Q. Similarly, 

Q_(t,r) = (l + (9(L / ))T_e +i ^ Lo(s)ds =W (t,r)- 1 + O(L / ). (16) 
We have for E(t) in the Born approximation: 

E(t)=i/ V (t,r)QL I (r)VU (t,r)- 1 dT. (17) 
Jo 

Finally, after making use of Eq. fll5ap . we get the ME for p in the Born approximation, 
which implies the neglect of the back-effect of the system on the bath, consistently with the 
weak-coupling assumption: 

p(t) = -i(h S (t) + (L I ) B )p(t) - f\h I {t)m Q {t,T)h I {T)VU Q {t,TY 1 ) B dT P {t) (18) 

Jo 

for Qp to t(0) = 0. 

C. Explicit equations for factorizable interaction Hamiltonians 

We now wish to write the ME explicitly for time-dependent Hamiltonians of the following 
form 12511 : 



H(t) = Hs(t) + H B + Hj(t), 
H T (t) = S(t)B; 



(19a) 
(19b) 



where Hs and H B are the system and bath Hamiltonians, respectively; and Hi, the inter- 
action Hamiltonian, is the product of operators S, B which act on the system (resp. bath) 
alone. We assume (Li)b = 0, by virtue of (B) B = 0. This also implies VLjV = 0. Equation 
( TT5j) now simplifies to: 



p(t) = -%L s {t)p{t) - [ rfr(L / (t)W (t,r)L / (r)PW (t,r)- 1 ) B p(t). 

Jo 



(20) 



Let us now write out the action of the operator Uq in terms of the unitary evolution 
operators of the system and bath: 



U {t, r)A = U s (t, r)U B (t - r)AU B (r - t)U s {t, r) t 



(21a) 



U s (t,r) = T + e-^rHs(t')dt' 



(21b) 



U B (t) = e~ lHBt . (21c) 
The integrand of the second term in Eq. (j20l . can now be written explicitly as: 



I(t,r) = Tr B [S(t)BMo(t,r) [S^BM^r)- 1 p{t) PB \] 



Tr 



B 



S(t)B, S{t,r)B{t-r),p{t)p B 



(22) 



where 5* and B are defined as: 



S(t,r) = Us(t,T)S(r)U s (t,T)\ 
B(t) = U b (t)BU b (-t). 



(23) 
(24) 



Using the commutativity of S and B, and of p B and H B , as well as the cyclic property of 
the trace, this gives after some rearrangement: 



I{t, r) = (BB(t - t)) b [S(t), S{t, r)p(t) 
Finally, defining the correlation function for the bath, 

$ T (t) = (BB(t)) B , 



+ H.c. 



(25) 



(26a) 



10 



we obtain the ME for p in the Born approximation: 



p(t)=-i[H s ,p(t)}+ J dr {<5> T (t - r) [s(t,r)p(t),S(t)\ + H.c.} . (26b) 

III. GENERALIZED BLOCH EQUATIONS 

Having derived the master equation, we focus on two regimes: a two-level system coupled 
to either an amplitude- or phase-noise (AN or PN) thermal bath. The bath Hamiltonian 
(in either regime) will be explicitly taken to consist of harmonic oscillators and be linearly 
coupled to the system (generalizations to other baths and couplings are obvious): 

H B = J2\ UJ x a x a \ ( 27 ) 
B = Zx(Kxa x + K* x a[). (28) 

Here a\, a\ are the annihilation and creation operators of mode A, respectively, and k x is 
the coupling amplitude to mode A. 

We use different modulation schemes for each regime, namely, dynamical off-resonant 
fields for the AN regime and time-dependent resonant fields for the PN regime. We derive 
the generalized Bloch equations for the two cases. 

A. Two-level system coupled to a thermal amplitude-noise bath 

We first consider the AN regime of a two-level system coupled to a thermal bath. We 
will use off-resonant dynamic modulations, resulting in AC-Stark shifts. The Hamiltonian 
(Eq. (|19a[) ) then assumes the following form: 

H s (t) = (u a + 5 a (t))\e)(e\ (29) 
S(t) = e(t)a x (30) 

where 5 a (t) is the dynamical AC-Stark shifts, e(t) is the time-dependent modulation of the 
interaction strength, and the Pauli matrix a x = \e)(g\ + \g)(e\. 

We derive the Bloch equations for the explicit case discussed above. Inserting Eqs. (12"9~|) - 
([30D into Eq. (l2lb|) and Eq. (J23]), we get: 

U s (t,T) = e -^--)-J>i<M*i)| e )( e | + ( 31 ) 
S(t,r) = e - lU}a(t - T) ' l ^ dhSa{tl) e{r)\e)(g\ + H.c. (32) 

11 



Plugging this into the ME (l26bl) . we arrive at the following modified Bloch equations: 

Pee = -Pgg = ~ Re{t) Pee + Rg{t) Pgg (33) 
Peg = Pge = ~ {(#(*) + lA (t)) + l[u a + 8 a (t)]} Peg 

+ {R(t) - lA a (t)} Pge, (34) 

where 

R(t) = [R e (t) + R g (t)}/2 (35) 

A a (t) = A e (t) - A g (t) (36) 

R e{g) (t)/2 + iA e(g) (t) = f dt'^t-t^K^t'y^-^ (37) 

lfe(*,0 = lC;(« ) = c(*K(0 (38) 
e(t) = i(t)e i ti dtlSa{tl) . (39) 

-R e (g) (^) is the modified downward (upward) transition rate of the excited (ground) state to 
the ground (excited) state. Their half-rate contributes to the decoherence rate, and A a (t) 
is the resonance (transition frequency) shift in energy due to the modified coupling to the 
bath. 

B. Two-level system coupled to thermal phase-noise bath 

Next, we consider the PN regime of a two- level system coupled to thermal bath, where we 
will use near-resonant fields with time- varying amplitude as our control. The Hamiltonians 
(Eq. p9aD ) then assume the following forms: 

H s (t) =u a \e)(e\ + V(t)a x (40) 
S(t) = e(t)a z (41) 

where V(t) = Vo(t)e~ luJat + c.c is the time-dependent resonant field, with real envelope Vo(t), 
e(t) is the time-dependent modulation of the interaction strength, a z = \e)(e\ — \g){g\- 

Since we are interested in dephasing, phases due to the (unperturbed) energy difference 
between the levels are immaterial. We eliminate this dependence by moving to the rotating 
frame. To avoid the need to time-order the propagator of the system Hamiltonian we tilt 
the rotating frame to the time-dependent basis: 

I T) = 4 (e~ W |e) + \g)) | |) = -L ( e -**t| e ) _ |^)) (42) 

12 



In this frame, the system and bath Hamiltonians become: 

Hs(t) = ^a z (43) 
S(t) = e(t)a x (44) 

where" denotes the rotated and tilted frame, a z = | |)(T l — I I an d & x = | |)(| | + | |)(T I- 
We can now derive the Bloch equations for the PN regime discussed above, and demon- 
strate their analogy to their AN counterparts (|33l . (!34|) . To this end we insert Eqs. (1431- (141 
into Eq. (|21b]) and Eq. ([23]), to get: 



U s (t,r) = e-^ 1 ^)/ 2 ! T)(T | +e l />i^o(ti)/2| ^ | (45) 
S(t,r) = e- l -fr dt ^)i(t)\ T)U | + H.c. (46) 

Plugging this into the ME (I26bl) . we arrive at the following modified Bloch equations: 

ptt = -hi = - R i(t)Prt + R i( t )Pu ( 47 ) 

PU = ph = -{(R(t)+iA a (t))+iVa(t)/2}p n + {R(t)-iA a (t)} Plh (48) 

where 

R(t) = [R^t) + R l (t)}/2 (49) 

A a (t) = A T (t)-A ; (t) (50) 

R m (t)/2 + iA m (t) = J dt'Mt-t')K ni) (t,t') (51) 

K ] it,t') = Kl(t,t')=e(t)e*(t') (52) 

e(t) = i(t)e l ti dtlVo{tl \ (53) 

As can be clearly seen, these modified Bloch equations are completely analogous to their 
AN counterparts, Eqs. (!33l) . (l34|) . provided we change the basis as follows: 

e ^T 9^1- (54) 

Despite their analogy, Eqs. (|37|) and (!5"T!) are not identical, due to the use of the rotat- 
ing frame in the PN case. Nevertheless, this analogy underscores the universality of our 
approach. 
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C. Spectral domain representation 



For both AN and PN regime, one can have a more insightful representation of the modified 
rates by transforming them to the frequency domain. For the long time limits (see Sec. IVT) . 
one arrives at the form 

/oo 
duG T {±{u a + u))F t {uj) (55) 
-oo 

/oo 
dujG T (±uj)F t (uj) (56) 
-oo 

where the difference is due to the fact that we used the rotating and tilted frame in the PN 
regime. Here Gt{oj) is the temperature-dependent bath coupling spectrum given by 

/oo 
dt$ T (t)e iuJt . (57) 
-oo 

Introducing the control-field fluence Q(t), the spectral modulation F t (u) can be normalized 
to unity: 



Q{t) = Hdr\e{r)\\ (58) 
Q(t) 

where 



*tf") = ^ (59) 



<*(«0 = 4= fdre{r)e^ (60) 

V £K JO 

is the finite-time Fourier transform of e(t). 

One can consider a more specific scenario, namely, coupling to zero-temperature (T = 0) 
AN bath. The effects of the bath then amount to the decay of the excited state's population, 
which can be written as: 

P e (t) = exp[-R e (t)Q(t)}. (61) 

IV. MODULATION ARSENAL FOR AN AND PN 

Any modulation with quasi-discrete, finite spectrum is deemed quasiperiodic, implying 
that it can be expanded as 

e(t) = J2^e- iUkt (62) 
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where Vk (k = 0,±1, ...) are arbitrary discrete frequencies such that 

Wk — v k ,\ > Q V/c 7^ k', (63) 

where Q is the minimal spectral interval. 

One can define the long-time limit of the quasi-periodic modulation, when 

ttt > 1 and t > t c , (64) 

where t c is the bath-memory (correlation) time, defined as the inverse of the largest spectral 
interval over which G t {uj) and Gt{— oj) change appreciably near the relevant frequencies 
u a + v k . In this limit, the fluence is given by 

Q(t)^e c t e c = ^| efe | 2 , (65) 

k 

resulting in the average decay rate: 

i? e = 27r^|A fe | 2 GK + z/ fc ), (66) 
k 

A fc = e fe /e c . (67) 



A. Phase modulation (PM) of the coupling 

1. Monochromatic perturbation 
Let 

e (t) = e e- tAt . (68) 

Then 

R e = 2irG T (cu a + A), (69) 

where A = const, is a frequency shift, induced by the AC Stark effect (in the case of atoms) 
or by the Zeeman effect (in the case of spins). In principle, such a shift may drastically 
enhance or suppress R relative to the Golden - Rule decay rate, i.e. the decay rate without 
any perturbation 

R GR = 2nG T (cu a ). (70) 

Equation 3.2 provides the maximal change of R achievable by an external perturbation, since 
it does not involve any averaging (smoothing) of G(uS) incurred by the width of F t (u>): the 

15 



modified R can even vanish, if the shifted frequency uj a + A is beyond the cutoff frequency 



of the coupling, where G(uj) = 



decoupling [10 



11. 



12 



13 
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Figure 
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i). This would accomplish the goal of dynamical 



181 ]. Conversely, the increase of R due to a shift 



can be much grea ter than that achievable by repeated measurements, i.e. the anti-Zeno 



effect 
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491 ] . In practice, however, AC Stark shifts are usually small for (cw) 



monochromatic perturbations, whence pulsed perturbations should often be used, resulting 
in multiple i>k shifts, as per Eq. fl66l) . 



2. Impulsive phase modulation 

Let the phase of the modulation function periodically jump by an amount <fi at times 
r, 2r, .... Such modulation can be achieved by a train of identical, equidistant, narrow 
pulses of nonresonant radiation, which produce pulsed AC Stark shifts of u a . Now 

e(t) = e i[i/T] *, (71) 

where [. . .] is the integer part. One then obtains that 



Q{t)=t, e c = l, (72) 

2 sin 2 (a;r/2) sin 2 [n(0 + lot) /2] 
nnruj 2 sin 2 [(0 + ur)/2] 

The excited-state decay, according to equation 55, has then the form (at t = nr) 



z. bin yuji i bin [ityip t wi . 



P e (nr) = exp[-R e (nr)nr], (74) 



where R e (nr) is defined by Eqs. 53 and 87. 

For sufficiently long times (Eq. fl64|) ) one can use Eq. fl66l) . with 



2/cvr <f> 4 sin 2 (0/2) 

Pk = — ~P |Afe| = (2kn - <Pf (?5) 



For small phase shifts, <f) <C 1, the k = peak dominates, 



whereas 



|Ao| 2 -l-g, (76) 



l^-^L (MO). (77) 



47r 2fc2 
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In this case one can retain only the k = term in Eq. (166|) . unless G{oS) is changing very 
fast with frequency. Then the modulation acts as a constant shift, (Fig. [T^,) 



A 



'T. 



(78) 



As |0| increases, the difference between the k = and k = 1 peak heights diminishes, 
vanishing for = ±7r. Then 



|A r 



|Ax| 2 = 4/tt 2 



(79) 



i.e., F t (u) for = ±7r contains two identical peaks symmetrically shifted in opposite directions 
(Figure[Tb) [the other peaks |Afc| 2 decrease with k as (2k — 1)~ 2 , totaling 0.19]. 

The foregoing features allow one to adjust the modulation parameters for a given scenario 
to obtain an optimal decrease or increase of R. Thus, the phase-modulation (PM) scheme 
with a small is preferable near a continuum edge (Figure Hk,b), since it yields a spectral 
shift in the required direction (positive or negative). The adverse effect of k ^ peaks in 
F t (uj) then scales as 2 and hence can be significantly reduced by decreasing |0|. On the 
other hand, if uj n is near a symmetric peak of G(u), R is reduced more effectively for ~ tt, 



as in Refs. 



11 



321 ]. since the main peaks of F t (uj) at ujq and U\ then shift stronger with r 1 
than the peak at ojq = —(p/r for 1. 

(b) 




jkF t (co) 



(C) 

G(«yT\ 



ft 



FIG. 1: Spectral representation of the bath coupling, G(lu), and the modulation, F t (u). (a) 
Monochromatic modulation, or impulsive phase modulation, with small phase shifts, <C 1, and 
1/r repetition rate, (b) Impulsive phase modulation, (-zr-pulses), = tt. (c) On-off modulation, 
with 1/ti repetition rate for t\ <C tq. 
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B. Amplitude modulation (AM) of the coupling 



Amplitude modulation (AM) of the coupling may be applicable to certain AN or PN 
scenarios. It arises, e.g., for radiative-decay modulation due to atomic motion through a 
high-Q cavity or a photonic crystal jsol, I5I] or for atomic tunneling in optical lattices with 



time- varying lattice acceleration [29, 



52|. 



1. On-off modulation 

The simplest form of AM is to let the coupling be turned on and off periodically, for the 
time 7~i and r — n, respectively, i.e., 

e(i ) = ( 1 i0InT ' < * < nT " + Tl ' (80) 
I for nro + n < t < (n + 1)tq 

(n = 0, 1, . . . ). Now Q(t) in 55 is the total time during which the coupling is switched on, 
whereas 

2sin 2 (u;ri/2) sin 2 (^r /2) 

i'nroK^) = 2 ■ 2/ 7^ ' 

7mrio; z sm (a;ro/2J 

so that 

P(nr ) =exp[-i2(7ir )nri], (82) 

where R(nro) is given by Eqs. 53 and 91. This case is also covered by 75, where the 
parameters are now found to be 

2 T l 2 7~l 2 f k-KTi 



e c = — , ^fc = , |A fc | = —sine . (83) 

ro r r V r o / 

It is instructive to consider the limit wherein t\ r and r is much greater than the 
correlation time of the continuum, i.e., G[iS) does not change significantly over the spectral 
intervals (27ik/r , 2tt (k + 1) / t ) . In this case one can approximate the sum 75 by the integral 
53 with 

F t (u) w (r 1 /27r)sinc 2 (^r 1 /2), (84) 
characterized by the spectral broadening ~ 1/ti (figure [Tb). Then equation 53 for R reduces 



23) 



to that obtained when ideal projective measurements are performed at intervals T\ 
Thus the AM on-off coupling scheme can imitate measurement-induced (dephasing) effects 
on quantum dynamics, if the interruption intervals ro exceed the correlation time of the 
continuum. 
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V. MULTIPARTITE DECOHERENCE CONTROL 



Multipartite decoherence control, for many qubits coupled to thermal baths, is a much 
more challenging task than single-qubit control since: (i) entanglement between the qubits is 
typically more vulnerable and more rapidly destroyed by the environment than single qubit 
coherence (ii) the possibility of cross- decoherence, whereby qubits are coupled 

to each other through the baths, considerably complicates the control. We have recently 



analyzed this situation and extended [30j, |37| the decoherence control approach of Sec. HIHlVl 



to multipartite scenarios, where the qubits are either coupled to zero-temperature baths or 
undergoing proper dephasing. 



A. Multipartite AN control by off-resonant modulation: singly excited systems 
coupled to T = baths 

The decay of a singly excited multi-qubit system (under amplitude noise) to the ground 
state, in the presence of off-resonant modulating fields is described by the following relaxation 



matrix 



30 



33| : 



oc 



J jf (t) = 2tt / duG jf (u)F tdj ,(u) (85) 

J — oo 

Gjji(u) = fir 2 ^ Vk,jVi,f8{u - w*) (86) 
k 

F t,jj>(^) = ~ w iKi'(^ - u r ) (87) 

Here Gjji(u) is the coupling spectrum matrix given by nature and F t .jj>(u>) is the dynamical 
modulation matrix, which we design at will to suppress the decoherence. The diagonal 
elements of the decoherence matrix are the time-integrated individual qubits' decay rates, 
while the off-diagonal elements are the cross-relaxation rates, pertaining to the coupling of 
the different qubits through the bath: virtual emission into the bath by qubit j and its 
virtual reabsorption by qubit j'. 

As an example, we may control the relaxation matrix elements by local (qubit-addressing) 
impulsive phase modulation, (see Sec. IIV A~2|) . described by 

iM = At ; ; -■ (88) 
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Here [...] denote the integer part, Tj and <fij are the pulse duration and the phase change 
for particle j, respectively. In the limit of weak pulses, of area <C n, Eq. (1BB1 yields 
e t j(uj) = e t j5{oj — Aj), where Aj = 4>jl T j is the effective spectral shift caused by the pulses. 

One can define the fidelity, F(t), total excitation probability, F p (t), and the autocorrela- 
tion function, F c (t) as follows: 

F(t) = Tr {]] (p(O)p(t)) (89) 

F p (t)=Tr {j ,#y( j (e\p(t)\ e h) ( 90 ) 
F c (t) = F(t)/F p (t) (91) 

where Tr^-j denotes tracing over all qubits. In the absence of dynamical control, the au- 
tocorrelation decays much faster than the total excitation probability, and is much more 
sensitive to the asymmetry between local particle-bath couplings. 

Thus, for initial Bell singlet and triplet states, which do not experience cross-decoherence 
but only different local decoherence rates, we find: 

\*(0)) = l/V2(\g) A \e) B ±\e) A \g) B ), (92) 

F p (t) = (e~ 2J ^ + e- 2J *«)/2; (93) 

F c (t) = (1 + C(t))/2 = 1/2 + e" AJ W/(l + e" 2AJ W), (94) 

AJ{t) = J A (t) - J B (t) (95) 



where C(t) is the concurrence 531 ] . 



Without any modulations, decoherence in this scenario has no inherent symmetry. Our 
point is that one can symmetrize the decoherence by appropriate modulations. The key is 
that different, "local", phase-locked modulations applied to the individual particles, accord- 
ing to Eq. (IH71) . can be chosen to cause controlled interference and/or spectral shifts between 
the particles' couplings to the bath. The F t jji{uj) matrices (cf . (j87]) ) can then satisfy 2N re- 
quirements at all times and be tailored to impose the advantageous symmetries described 
below. By contrast, a "global" (identical) modulation, characterized by F t> jj>(uj) = |et(a;)| 2 , 
is not guaranteed to satisfy N ^> 1 symmetrizing requirements at all times (Fig. [2^). 

The most desirable symmetry is that of identically coupled particles (ICP), which would 
emerge if all the modulated particles could acquire the same dynamically modified deco- 
herence and cross-decoherence yielding the following N x N fully symmetrized decoherence 
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(a) 

\e tJ (a>)?=\e,j.(a>)f 



(b) 




FIG. 2: Two two-level particles in a cavity, coupled to the cavity modes (thin lines) and subject to 
local control fields (thick lines). (a,b) Frequency domain overlap of coupling spectrum (dotted) and 
modulation matrix elements (solid), resulting in modified decoherence matrix elements (shaded), 
for: (a) global modulation (ICP symmetry), (b) cross-decoherence elimination (IIP symmetry), (c) 
General modulation scheme. 

matrix 

Jjf(t)=r(t) Vj,f. (96) 

ICP would then give rise to a (N — l)-dimensional decoherence-free subspace: the entire 
single-excitation sector less the totally symmetric entangled state. An initial state in this 
DFS [19] would neither lose its population nor its initial correlations (or entanglement). 

Unfortunately, it is generally impossible to ensure this symmetry, since it amounts to 
satisfying iV(iV — 1)/2 conditions using N modulating fields. Even if we accidentally succeed 
with iV particles, the success is not scalable to iV + 1 or more particles. Moreover, the ability 
to impose the ICP symmetry by local modulation fails completely if not all particles are 
coupled to all other particles through the bath, i.e. if some Gjj/(uj) elements vanish. 

A more limited symmetry that we may ensure for N qubits is that of independent identical 
particles (IIP). This symmetry is formed when spectral shifts and/or interferences imposed 
by iV modulations cause the N different qubits to acquire the same single-qubit decoherence 
r(t) and experience no cross-decoherence. To this end, we may choose £t,j(u) ~ e t jS(u — Aj). 
We shall deal with N identical qubits, and set Uj = ojq. We also require that at any chosen 
time t — T, the AC Stark shifts satisfy J Q T rfr5j(r) = 2irm, where m = 0, ±1, .... This 
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requirement ensures that modulations only affect the decoherence matrix f|85l) . but do not 
change the relative phases of the entangled qubits when their MES is probed or manipulated 
by logic operations at t = T. 

The spectral shifts Aj can be different enough to couple each particle to a different 
spectral range of bath modes so that their cross-coupling vanishes: 

Jjf(t) = 2ne* J e tJ/ J dujG j:j >(ujQ + oj)5(uj - Aj)5(uj - A f ) -> 0. (97) 

Here, the vanishing of Gjj>(u) for some j,j' is not a limitation. The N single-particle 
decoherence rates can be equated by an appropriate choice of TV parameters {A,,}: 

Jgf(t) = 27r|e t , i | 2 G i > + A,) = 5 jf r(t), (98) 

where Sjf is Kronecker's delta (Fig. [Tb). The IIP symmetry results in complete correlation 
preservation, i.e. F c (t) = 1, but still permits excited-state population loss, F p (t) = e ~ 2Re 'f r '(*^ 
(Fig. [3]). If the single-particle r(t) may be dynamically suppressed, i.e. if the spectrally 
shifted bath response Gjj(uOj + Aj) is small enough, this F p (t) will be kept close to 1. 




t 

FIG. 3: Fidelity of the IIP symmetry for two TLS coupled to zero-temperature baths. The initial 
state is entalged, |^(0)) = l/\/2 {\g) A \e) B + \e) A \g) B ). 

B. Multipartite PN control by resonant modulation 

One can describe phase-noise, or proper dephasing, by a stochastic fluctuation of the 
excited-state energy, uj a — > uj a + S r (t), where S r (t) is a stochastic variable with zero mean, 
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(8(t)) = 0, and (5(t)8(t')) = $ p (t — t') is the second moment. For multipartite systems, 
where each qubit can undergo different proper dephasing, 5 jit), one has an additional sec- 
ond moment for the cross-dephasing, (8j(t)8f(t')) = §fj/(t — t'). A general treatment of 



30| . Here we 



multipartite systems undergoing this type of proper dephasing is given in Ref. 
give the main results for the case of two qubits. 

Let us take two TLS, or qubits, which are initially prepared in a Bell state. We wish to 
obtain the conditions that will preserve it. In order to do that, we change to the Bell basis, 
which is given by 

\B l , 2 ) = l/V2e^ t {\e) l \g) 2 ±\g) l \e) 2 ) (99) 
\B 3A ) = 1/V2 (e^|e)i|e> 2 ± {g)^),) . (100) 

For an initial Bell-state p z (0) = \Bi){Bi\, where / = 1...4, one can then obtain the fidelity, 
F l {t) = {B l \p l {t)\B l ), as: 

F l {t) = cos(0 ± (t))Re [e^W (l - \ £ .., Jf f> { (t))] , (101) 

where 

hit) = 2 £ drVoAr) (102) 

Jf ftl (t) = 2tt dLjG^FwM (103) 

G fA") = S-oo^fA^ ( 104 ) 

%(«) = Mw)| J (105) 

i-f.jrM = -F t , jr ,M = 4,M e iM ( 106 ) 

F tjj'A u ) = - F t,n>M = e t j{u)ei r {uj) (107) 

where V j{t) is the amplitude of the resonant field applied on qubit j, <f>±(t) = {<fii{t) ± 
2 (^))/2 and the + corresponds to k = 1,3 and 0_ to k = 2,4. Expressions fll 1 [) - (11 7j) 
provide our recipe for minimizing the Bell-state fidelity losses. They hold for any dephasing 
time-correlations and arbitrary modulation. 

One can choose between two modulation schemes, depending on our goals. When one 
wishes to preserve and initial quantum state, one can equate the modified dephasing and 
cross-dephasing rates of all qubits, Jjj>,i(t) = J if). This results in complete preservation of 
the singlet only, i.e. F 2 {t) = 1, for all t, but reduces the fidelity of the triplet state. On 
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the other hand, if one wishes to equate the fidelity for all initial states, one can eliminate 
the cross-dephasing terms, by applying different modulations to each qubit (Fig. 0]), causing 
F t jjij(uj) = Vj 7^ j'. This requirement can be important for quantum communication 
schemes. 




FIG. 4: Cross-decoherence as a function of local modulation. Here two qubits are modulated by 
continuous resonant fields, with amplitudes f2i,2- The cross-decoherence decays as the two qubits' 
modulations become increasingly different. The bath parameters are $t(*) = e _t// * c , where t c = 0.5 
is the correlation time; and fii = 3. 



VI. CONCLUSIONS 

In this paper we have expounded our universal approach to the dynamical control of 
qubits subject to AN and PN, by either off- or on-resonant modulating fields, respectively. 
It is based on a general non-Markovian master equation valid for weak system-bath coupling 
and arbitrary modulations, since it does not invoke the rotating wave approximation. The 
resulting universal convolution formulae provide intuitive clues as to the optimal tailoring of 
modulation and noise spectra. Our analysis of multiple, field-driven, qubits which are cou- 
pled to partly correlated or independent baths or undergo locally varying random dephasing 
has resulted in the universal formula (|85|) for coupling to zero-temperature bath, and (lIOip 
for Bell-state preservation under local proper dephasing. 

Our general analysis allows one to come up with an optimal choice between global and 
local control, based on the observation that the maximal suppression of decoherence is not 
necessarily the best one. Instead, we demand an optimal phase-relation between different, 
but synchronous local modulations of each particle. The merits of local vs. global modula- 
tions have been shown to be essentially twofold: 

• Local modulation can effectively decorrelate the different proper dephasings of the 
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multiple TLS, resulting in equal dephasing rates for all states. For two TLS, we 
have shown that the singlet and triplet Bell-states acquire the same dynamically- 
modified dephasing rate. This should be beneficial compared to the standard global 
"Bang-Bang" (7r-phase flips) if both states are used (intermittently) for information 
transmission or storage. 

• For different couplings to a zero-temperature bath, one can better preserve any initial 
state by using local modulation which can reduce the decay as well as the mixing with 
other states, than by using global modulation. It was shown that local modulation 
which eliminates the cross-decoherence terms, increases the fidelity more than the 
global modulation alternative. For two TLS, it was shown that local modulation 
better preserves an initial Bell-state, whether a singlet or a triplet, compared to global 
7r-phase "parity kicks" . 
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